Observation of an enhancement in e+e — t- T(15')7r+7r , T(2S')7r+7r , and T(3S')7r+7r 

production near y/s = 10.89 GeV 
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We measure the production cross sections for e^e~ — s> T(lS)7r+7r^, T(25')7r'''7r~ , and T(3S)7r+7r~ 
as a function of ^/s between 10.83 GeV and 11.02 GeV. The data consists of 8.1 fb~^ collected with 
the Belle detector at the KEKB e'^e" collider. We observe enhanced production in all three final 
states that does not agree well with the conventional T( 10860) lineshape. A fit using a Breit- 
Wigner resonance shape yields a peak mass of [10888.4^2.6 (stat) ± 1.2(syst)] MeV/c^ and a width 
of [30.7l7.o(stat) ±3.1(syst)] MeV/c^ 

PACS numbers; 13.25.Gv, 14.40.Pq 



A host of new charmonium-like mesons have been dis- 
covered recently that do not seem to fit into the con- 
ventional cc spectrum. Possible interpretations for these 
exotic states include multiquark states, cq-cq mesonic- 
molecules (where q represents a u-, d- or s-quark), or ccq 
"hybrids" (where g is a gluon). The narrow X(3872) 
and various vector mesons, in particular, the broad 
y(4260) [H3 resonance, were revealed by their dip- 
ion transitions to J/ip (and ip'). By analogy with the 
y (4260) discovery, it was suggested that a hidden-beauty 
counterpart, denoted Yb Q, could be sought in radiative 
return events from e'^e~ at center-of-mass energy on the 
T(10860) peak, or by an energy scan above it. 

Analyzing a data sample recorded by the Belle detec- 
tor at a single energy near the T(10860) resonance peak, 
no unusual structure was observed below the peak. In- 
stead, anomalously large T(15)7r+7r~ and T(2S')7r+7r~ 
production rates were observed [a] . If these signals are at- 
tributed entirely to dipion transitions from the T(10860) 
resonance, the corresponding partial widths would be 
two orders of magnitude larger than those for corre- 
sponding transitions from the T{2S), T(35), and T(4S') 
states. Possible explanations include a new nonperturba- 
tive approach [7| to the decay widths of dipion transitions 
of heavy quarkonia, the presence of final state interac- 
tions ^] , or the existence of a tetraquark state [1, [l^l ■ A 
measurement of the energy dependence of the cross sec- 
tions for e+e~ — T(nS')7r+7r~ (n = 1, 2, 3) in and above 
the T(10860) energy region provides more information 
for exploring these models. 

Here we report the observation of the production of 
e+e- -> T(15)7r+7r-, T(25)7r+7r-, and T(35)7r+7r- at 
~ 10.83, 10.87, 10.88, 10.90, 10.93, 10.96, and 11.02 
GeV based on a previous data sample of 21.7 fb~^ already 
reported in Ref. ^] and a new 8.1 fb~^ data sample at 
the other energies, both collected with the Belle detector 



at the KEKB e+e collider [ll|. In addition, nine scan 



points with an integrated luminosity of ~ 30 pb^^ each, 
collected in the range of y^s = 10.80-11.02 GeV, are used 
to obtain the hadronic line shape of the T(10860). 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector, a cen- 
tral drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel-like arrangement of 
time-of-flight scintillation counters, and an electromag- 
netic calorimeter (ECL) comprised of CsI(Tl) crystals 
located inside a superconducting solenoid that provides 
a 1.5 T magnetic field. An iron flux- return located out- 
side the coil is instrumented to detect mesons and 
to identify muons (KLM). The detector is described in 
detail elsewhere (l^ . 

Events with four well-reconstructed charged tracks and 
zero net charge are selected to study T(nS')7r+7r~ pro- 
duction. A final state that is consistent with a T{nS) 
candidate and two charged pions is reconstructed. A 
T(nS) candidate is formed from two muons with op- 
posite charge, where the muon candidates are required 
to have associated hits in the KLM detector that agree 
with the extrapolated trajectory of a charged track pro- 
vided by the drift chamber. The T{nS) — >■ e+e~ chan- 
nels have little sensitivity due to the lower dielectron 
efficiency and the 20 times larger background from ra- 
diative Bhabha events. The other two charged tracks 
must have a low likelihood of being electrons (the like- 
lihood is calculated based on the ratio of ECL shower 
energy to the track momentum, dE/dx from the CDC, 
and the ACC response); they are then treated as pion 
candidates. This requirement suppresses the background 
from e+e^ — >■ /i+/i~7 — ^ /i+/i~e+e~ with a photon con- 
version. The cosine of the opening angle between the tt"*" 
and TT^ momenta in the laboratory frame is required to 
be less than 0.95. The four-track invariant mass must 
satisfy |M(/i+//-7r+7r-) - v^l < 150 MeV/c^. The trig- 
ger efficiency for four-track events satisfying these criteria 
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TABLE I: Center-of-mass energy (-y/s), integrated luminosity [C), signal yield (A^s), reconstruction efficiency, and measured 
cross section (cr) for e^e~ — >■ T(lS')7r"'"7r~, T(2S)7r'*'7r~, and T(3S')7r^7r~ . Due to a negative yield, an upper limit at 90% 
confidence level for cr[e+e~ — ^ T(3S')7r"'"7r^] at = 10.9555 GeV is given. 
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(a)Y(1 S)3t7t,\ s=1 0.8255GeV 



(b)Y(1 S)7tJt,\ s=1 0.8805GeV 



(c)Y(1S)itJt,\s=10.8955GeV 
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L (g)Y(2S)it7t,\ s=10.8255GeV 




(h)Y(2S)ii7i,\s=10.8805GeV 



(i)Y(2S)it7t,\'s=1 0.8955GeV 
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(k)Y(2S)7t7t,\ s=1 0.9555GeV 
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FIG. 1: The distributions of AA/~[Vi-MT{„s)] (n = 1,2,3) for (a-f) T(lS')7r+7r- , (g-1) T(2S)7r+7r-, and (m-r) T(35')7r+7r- 
events with the fit results superimposed. The six columns of plots represent the data samples collected at different energies. 
The dashed curves show the background components in the fits. 



is very close to 100%. 

The kinematic variable AM, defined by the differ- 
ence between M(/Lt+/z~7r+7r~) and M(/z"'"/Lt~), is used 
to identify the signal candidates. Sharp signal peaks 



will occur at AA/ = — Afx(nS)- The candidate 
events are separated into three distinct regions defined 
by I AM - [Vs - Mt(„s)]| < 150 MeV/c^ for n = 1, 2, 
and 3, and signal yields are extracted from an unbinned 



4 



extended maximum likelihood fit to the AM distribution 
within each region. The likelihood function for each fit 
is defined as 



-{N,+N,) N 

L{Ns,Nt) = Y[[N,-PsiAM,)+Nb-Pb{AM,)] , 

i=i 



where Ng (Nt) denotes the yield for signal (background), 
and Pg (Pb) is the signal (background) probability den- 
sity function (PDF). The signal is modelled by a sum of 
two Gaussians while the background is approximated by 
a linear function. The Gaussians parameterized for the 
signal PDF are fixed from the Monte Carlo (MC) simu- 
lation at each energy point. We fit 18 AM distributions 
simultaneously with common corrections to the mean and 
width of the signal Gaussians. The widths are found to 
be around 19 ±8% higher than the expectations given by 
the MC simulations. Figure [T] shows these AM distribu- 
tions with fit results superimposed; the other three AAI 
distributions for the data collected at ^/s « 10.867 GeV 
are included in Ref. @. 

The measured signal yields, reconstruction efficiencies, 
integrated luminosity, and the production cross sections, 
as well as the results from the previous publication Q 
for the data sample collected at ^/s = 10.867 GeV, are 
summarized in Table H] The efficiencies for T(3S')7r+7r^ 
are much improved compared to Ref. Q as the inefficient 
selection criterion ^max < 175° has been removed, where 
^max is the maximum opening angle between any pair of 
charged tracks in the center-of-mass frame. This channel 
is clearly seen for the first time. The selection criteria 
already constrain the reconstructed total energy to be 
very close to the input beam energy. The contribution 
from radiative return events is expected to be very low. 
The observed peaks in AM distributions agree with MC 
expectations without initial state radiation. 

For the cross section measurements, systematic uncer- 
tainties are dominated by the T{nS) — ^ branch- 
ing fractions, reconstruction efficiencies, and PDF pa- 
rameterization for the fits. Uncertainties of 2.0%, 8.8%, 
and 9.6% for the T(15'), T(25), and T{3S) ^i+fi- 
branching fractions are included, respectively. For the 
T(lS')7r+7r~ and T(2S')7r"'"7r~ modes, the reconstruction 
efficiencies are obtained from MC simulations using the 
observed M(7r+7r~) and cos6'ho1 (the angle between the 
7r~ and T(10860) momenta in the tt+tt^ rest frame) 
distributions in our previous publication as inputs @. 
The uncertainties associated with these distributions give 
rise to 2.7%-4.5% and 1.9%-4.2% uncertainties for the 
T(lS')7r+7r~ and T(2S')7r+7r~ efficiencies, respectively. 
The ranges on the uncertainty arise from the energy de- 
pendence of the 7r+7r~ system. We use the model of 
Ref. as well as a phase space model as inputs for 
T(3iS')7r"'"7r~ measurements; the differences in acceptance 
are included as systematic uncertainties. The uncertain- 
ties from the PDF parameterization are estimated either 
by replacing the signal PDF with a sum of three Gaus- 
sians, or by replacing the background PDF with a second- 



order polynomial. The differences between these alterna- 
tive fits and the nominal results are taken as the system- 
atic uncertainties. Other uncertainties include: track- 
ing efficiency (1% per charged track), muon identifica- 
tion (0.5% per muon candidate), electron rejection for 
the charged pions (0.1-0.2% per pion), trigger efficien- 
cies (0.1-5.2%), and integrated luminosity (1.4%). The 
uncertainties from all sources are added in quadrature. 
The total systematic uncertainties are 7%-ll%, 11%- 
16%, and 12%-14% for the T(lS')7r+7r-, T(2S')7r+7r-, 
and T(3S')7r+7r^ channels, respectively. 

In order to extract the resonance shape, we perform a 
fit to the measured production cross sections, includ- 
ing the one estimated at y/s = 10.867 GeV, using the 
model 



<^Arrrn\Ro + e"'>BW{fi,r)\ 

where crj]^ = Ana^ /3s is the leading-order e^e 
cross section and BW{ii,T) is the Breit-Wigner function 
l/[(s — /i^) + i/iF]. The normalizations for T(lS')7r+7r~, 
T(2S')7r+7r^, and T(3S')7r+7r~ (Axtttt), as well as the am- 
plitude of the fiat component i?o, the mean /i, width F, 
and the complex phase of the parent resonance are free 
parameters in the fit. Because of the low statistics, com- 
mon shape parameters (/i, F, (j), and i?o) are introduced 
for the three different final states. Results of the fits, 
shown as the smooth curves in Fig. [21 are summarized 
in Table HIl The fit quality is = 24.6 for 14 degrees 
of freedom (corresponding to a confidence level of 3.9%). 
An alternative fit without the last data point collected at 
^ 11.02 GeV yields a similar result, ^ = 10889.0t^J 



MeV/c2, F 



37t]n MeV/c2, and = 21.3 for 11 de- 



-10 



grees of freedom (corresponding to a confidence level of 
3.0%). Systematic uncertainties associated with the cross 
section measurements are propagated to the resonance 
shapes. The fits are repeated, and the variations on the 
shape parameters are included as the systematic uncer- 
tainties. In addition to the uncertainties on the cross 
sections, the beam energy around the T(10860) is mea- 
sured by Mx(nS) + in the T(n5')7r+7r~ events, and 
an uncertainty of ±1 MeV (comprising the uncertainties 
of Mx{nS) given in Ref. flTj and of AM in Ref. f^]) is in- 
cluded. For the scan data, a common energy shift is also 
obtained from the fit to T(n5')7r"'"7r~ events. The relative 
beam energies are further checked using the M{^'^fi~) 
distributions of /i-pair samples. 

We also determine the resonance parameters for the 
T(10860) using bb hadronic events from the energy scan 
at y/s between 10.80 and 11.02 GeV. We measure the 
fraction Rb = CTb/cr^^, where ab — N^^"^^'"^ {s) / Ceb{s) is 
the e+e^ — ^ bb hadronic cross section. The number of 
e+e" bb events with i?2 < 0.2 (7Vjf^<°-^) is estimated 
by subtraction of non-56 events scaled from a data set 
collected at ^/s ~ 10.52 GeV, where i?2 denotes the ra- 
tio of the second to zeroth Fox- Wolfram moments 
Selection criteria for hadronic events are described in 
Ref. 15|. The acceptance for e+e~ — bb (eb(s)) is found 
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FIG. 2: The energy-dependent cross sections for e e~ — )■ 
T(n5')7r+7r~ (n = 1, 2, 3) processes normalized to the leading- 
order e'^e~ — >■ n'^ fj,~ cross sections. The results of the fits are 
shown as smooth curves. The vertical dashed line indicates 
the energy at which the hadronic cross section is maximal. 



to vary slightly from 68.1% to 70.5% over the range of 
scan energies. The line shape used to model our data is 
given by \Anr\^ + \Ao+Aw86oe"'"°'"'° BW {nwseo, rio86o) + 
^ii02oe''^"°"'Siy(/xiio2o,riio2o)P; this parameteriza- 
tion is the same as that used in Ref. [l^. We perform a 
fit to our Ri, measurements as shown in Fig.[31[a). The 
shapes for T(11020) ((^nn20, M11020 , and riio2o) are fixed 
to the values in Ref. since our data points are not 
able to constrain the T(11020) parameters. The resulting 
shape parameters for the T(10860) are ^loseo = 2.33jlQ'24 
rad, ^10860 10879 ± 3 MeV/c^ , and Tioseo = 46i? 
MeV/c^. These values are consistent with those obtained 
in Ref. [H]. The quality of the fit is = 4.4 for 9 de- 
grees of freedom (corresponding to a confidence level of 



Figure [3jb) shows the ratio between cr[e+e^ ^■ 
T(nS')7r+7r^] and a[e~^e^ bb] as a function of y/s. 
As shown in Fig. [3{c), an alternative fit applied to the 
T(nS')7r+7r^ cross sections and Rb measurements, simul- 
taneously, with common values of /xttttt — /^loseo and 



TABLE II: Cross sections at peak (o-pcak), mean (/i), width 
(r), phase (0), and the amplitude for the constant com- 
ponent (Ro) from the fit to the energy-dependent e^e~ — >■ 
T(15)7r^7r~, T(25')7r^7r~ , and T(35')7r"'"7r~ cross sections. 
Normalization of the resonance term in the fitting model is 
represented by apeak. There are two solutions for <jf> and _Ro 
with identical x^- The first uncertainty is statistical, and the 
second is systematic. 
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^ peak 
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FIG. 3: (a) R^, and (b) the ratio between G\e^ e~ — >■ 
T(n5)7r^7r~] and (j\e^ e~ — >■ 66] as a function of ~/s\ 
(c) the energy-dependent cross section ratios for e" — > 
T(nS')7r^7r~ events, the result of fits with resonant param- 
eters from _Ri, or PDG averages are superimposed. The hor- 
izontal dotted line in (a) is the non- interfering |AnrP contri- 
bution in the fit. The vertical dashed line indicates the energy 
at which the hadronic cross section is maximal. 



Fttttt = Fioseo increases the by 8.71 and reduces the 
degrees of freedom by two. This corresponds to a de- 
viation of 2.5cr from the nominal fit with separated /i 
and r for T(n5')7r+7r~ and Rh measurements. We also 
examine the systematic sources as described above, and 
the smallest deviation obtained is 2.0(7. If the resonant 
mean and width are fixed to the results in Ref. [l^ or 
PDG values [13], deviations of 3.9cr or 5.6cr, respectively, 
are obtained. Scans within the region with (/^ttttt ~ 

M1O860) /o'(/il086o) ± (Fttttt — Fioseo) /O'(ri086o) < 1, 

yield minimum deviations of 3.4(t or S.lcr, respectively. 
In summary, we report the observation of enhanced 

+^-, 

10.83 and 11.02 GeV. The 



e+e- -> T(15)7r+7r-, T(2S')7r+7r-, and T(3S')7r+7r- pro- 
duction at i/s between \fs : 
energy- dependent cross sections for e" — >■ T(nS')7r~'"7r~ 
events are measured for the first time, and are found to 
differ from the shape of the €^e~ — >■ bb cross section. 
A Breit-Wigner resonance shape fit yields a peak mass 
of 10888.4^2-7 (stat) ± 1.2 (syst) MeV/c^ and a width of 

30.71:^;^ (stat) ± 3.1 (syst) MeV/c^. A fit excluding the 
-^i ^ 11.02 GeV data point is consistent with the nomi- 
nal fit, indicating no strong contribution of T(n5')7r"'"7r~ 
events from T(11020). The T(10860) shape parame- 
ters obtained from our R\y hadronic cross section mea- 
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surements are consistent with the measurements from 
BaBar. The differences between the shape parameters 
from T(riS')7r+7r~ events and from our i?b measurements 
are ^ttttt - Mioseo = 9 ± 4 MeV/c^ and Tttvtt - Tioseo = 
-15j:i2 MeV/c^. A fit to the e+e" ^ T(nS')7r+7r- 
cross sections with our measured T(10860) mean and 
width yields a statistical deviation of 2.5a (2. Oct includ- 
ing systematic uncertainties). The T(n5)7r+7r^ partial 
widths were found to be much larger than the expecta- 
tions for conventional T{5S) states. As an extension, 
energy- dependent T(nS')7r+7r^ production cross sections 
are measured; the observed structure deviates slightly 
from the T(10860) shape obtained from the hadronic 



cross sections. More data is required to clarify the reso- 
nance structure. 
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